QUANTIZING NON-LAGRANGIAN GAUGE THEORIES: 
AN AUGMENTATION METHOD 



S.L. LYAKHOVICH AND A.A. SHARAPOV 

Abstract. We discuss a recently proposed method of quantizing general non-Lagrangian 
gauge theories. The method can be implemented in many different ways, in particular, it 
can employ a conversion procedure that turns an original non-Lagrangian field theory in d 
dimensions into an equivalent Lagrangian, topological field theory in d+ 1 dimensions. The 
method involves, besides the classical equations of motion, one more geometric ingredient 
called the Lagrange anchor. Different Lagrange anchors result in different quantizations 
of one and the same classical theory. Given the classical equations of motion and La- 
grange anchor as input data, a new procedure, called the augmentation, is proposed to 
quantize non-Lagrangian dynamics. Within the augmentation procedure, the originally 
non-Lagrangian theory is absorbed by a wider Lagrangian theory on the same space-time 
manifold. The augmented theory is not generally equivalent to the original one as it has 
more physical degrees of freedom than the original theory. However, the extra degrees of 
freedom are factorized out in a certain regular way both at classical and quantum lev- 
els. The general techniques are exemplified by quantizing two non-Lagrangian models of 
physical interest. 



We are thankful to G. Barnich and M. Henneaux for useful discussions of this work, and we appreci- 
ate warm hospitality of International Solvay Institutes and Department of Mathematical and Theoretical 
Physics of University Libre de Bruxelles where this work has been finalized. The work was partially sup- 
ported by the RFBR grant 06-02-17352 and the grant for Support of Russian Scientific Schools 5103.2006.2 

and also benefited from the "Innovative Education Programme" of TSU. 

1 



2 



s.l. lyakhovich and a.a. sharapov 
Contents 



1. Introduction 3 

2. Lagrange structure and S'oo-algebras 7 

2.1. Classical dynamics 7 

2.2. Regularity conditions __|8 

2.3. Lagrange structure 

2.4. S'oo-algebras 

3. BRST complex 

3.1. An ambient symplectic supermanifold 

3.2. BRST charge 

3.3. Hamiltonian interpretation 

3.4. Physical observables 

4. Quantization 

4.1. Quantum BRST cohomology 

4.2. Generalized Schwinger-Dyson equations 

4.3. Path-integral representation 

5. Augmentation 

5.1. An augmented BRST complex. 

5.2. Interpretation 

5.3. Quantizing non-Lagrangian dynamics via augmentation 

6. Examples of quantizing non-Lagrangian field theories 

6.1. Maxwell electrodynamics with monopoles 

6.2. Self-dual p-form fields 
References 



1. Introduction 



3 



Classical dynamics can be consistently formulated in terms of equations of motion alone. 
The variational principle, being a useful tool for studying various aspects of classical dy- 
namics, is not needed to have the classical theory defined as such. However, to promote 
the classical dynamics to quantum level, it is insufficient to know only the equations of 
motion, one or another extra structure is needed. If the quantum theory is supposed to 
be formulated in the language of Feynman's path integrals, it is the action functional that 
can serve as the additional ingredient needed for quantization. On the other hand, any 
Lagrangian equations of motion can always be brought to a (constrained) Hamiltonian 
form that makes possible applying canonical quantization. Furthermore, the method of 
deformation quantization applies to the Hamiltonian systems even though the underlying 
Poisson bracket is degenerate [1] (in which case the Hamilton equations can have no vari- 
ational formulation). As it has been recently found [2], [3], the deformation quantization 
can also be implemented under far less restrictive conditions on the equations of motion 
than the requirement to be Hamiltonian. Roughly speaking, the phase-space evolution flow 
is not required to be Hamiltonian: It is sufficient if the evolution preserves the Poisson 
bracket modulo constraints and gauge transformations. The bracket, in its turn, is also 
required to satisfy the Jacobi identity in a weak sense, i.e., modulo constraints and gauge 
transformations. For accurate definitions, see [2], [3]. 

So, the deformation quantization has progressed in recent years reaching far beyond the 
range of theories admitting variational principle for equations of motion. At the same time, 
the methods of constructing the partition function^j dating back to Feynman, Schwinger 
and Dyson, and being now developed in full generality for arbitrary Lagrangian gauge 
theories |4], [5], [6], have not made much progress in the class of theories having no action 
functional. Until recently, no general method has been known to path-integral quantize a 
non-Lagrangian theory as it was not clear what might be a generalization of the familiar 
Schwinger-Dyson's equation in the situation where no Lagrangian formulation is possible 
for the classical equations of motion. 

In our recent papers [7j, [8], we have identified a general structure, called the Lagrange 
anchor, which is determinative for the quantization in terms of partition functions in the 
same sense as the Poisson bracket defines deformation quantization in terms of a star 
product. The Lagrange anchor is a geometric object that can be interpreted in many 
different ways. In particular, one could say that the Lagrange anchor is related to the 



^In the following we will also use the term probability amplitude as a synonym for the partition function. 
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canonical anti-bracket of the Batalin-Vilkovisky formalism [5] much like a generic (i.e., 
possibly degenerate and non-constant rank) Poisson bracket is related to the canonical 
Poisson bracket. The anchor is also required to satisfy certain compatibility conditions 
involving equations of motion. In Lagrangian theory, these conditions are automatically 
satisfied for the canonical anti-bracket in consequence of the fact that the equations of 
motion are variations of the action functional. If the anchor is invertible, these compatibility 
conditions ensure existence of the equivalent Lagrangian formulation. It turns out that the 
partition function can still be constructed by making use only of the equations of motion and 
the Lagrange anchor, even though the latter is degenerate, defining no action functional. 
The next section contains an accurate definition of the Lagrange anchor and discussion of 
its properties. 

By now two procedures have been worked out to construct partition functions for gen- 
eral non-Lagrangian theories. The first one [7] suggests a conversion of the original non- 
Lagrangian field theory in d dimensions into an equivalent {d + l)-dimensional Lagrangian 
topological theory which can then be quantized by the standard BV method. The conver- 
sion procedure is quite ambiguous and essentially depends on the choice of the Lagrange 
anchor. If the anchor was invertible (that assumes implicitly the existence of some action 
for the original dynamics), the path integral can be explicitly taken in the bulk of the 
topological theory resulting in Feynman's partition function for original action. With a 
general (non-invertible) anchor, the answer for the partition function cannot be reduced 
to the canonical form, but it remains fully consistent and allows quite natural physical 
interpretation [8]. If the anchor is chosen to be zero, the partition function will correspond 
to the classical transition amplitude [9J. The second method to quantize a classical the- 
ory with non-Lagrangian equations of motion [8] suggests a nontrivial generalization of the 
Schwinger-Dyson equation that any partition function must satisfy. This equation, involv- 
ing classical equations of motion and the Lagrange anchor, reduces to the BV quantum 
master equation whenever the anchor is invertible. 

In this paper, we propose an alternative procedure of constructing partition functions for 
general dynamical systems. This procedure starts with the same input data: classical equa- 
tions of motion and the Lagrange anchor, but it exploits quite different idea and technology. 
We call this procedure an augmentation because it is motivated by a widespread view that 
either a non-Lagrangian system can be reshaped into an equivalent Lagrangian model in an 
appropriately extended configuration space, or it describes an effective dynamics emerging 
from a Lagrangian theory after averaging over some degrees of freedom or their exclusion 
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from the equations of motion. So, the intuitive intention about quantizing a non-Lagrangian 
theory is to augment it first to a Lagrangian one, and then the augmented theory can be 
quantized in the usual way. No general method is known to date to equivalently reformu- 
late any given non-Lagrangian model as a Lagrangian one by adding a finite number of 
new fields. We propose a uniform procedure to construct an augmented Lagrangian the- 
ory for any (non-) Lagrangian dynamics, which is not however an equivalent reformulation. 
The augmented theory may have, in principle, more degrees of freedom than the original 
model, but classically, the original dynamics are easily singled out by imposing appropriate 
boundary conditions on the extra fields. These boundary conditions guarantee that the 
original fields evolve precisely in the same way as in the original theory, while new fields do 
not evolve at all. This reduction mechanism always restores the original dynamics in the 
augmented theory including the case where the original theory is Lagrangian. Quantizing 
the augmented Lagrangian system by conventional BV procedure and integrating the new 
fields out in the path integral, one gets the original (not necessarily Lagrangian) dynamics 
quantized. If the original theory is Lagrangian, the integral can be taken explicitly over the 
augmentation fields with corresponding boundary conditions, and the partition function 
obtained in this way will coincide with that constructed from the BV master action for 
the original Lagrangian. If the original theory is not Lagrangian, the constructed partition 
function is still correct that can be seen in several ways, although it cannot be represented 
anymore as an exponential of any (local or non-local) action functional. 

Let us also comment on an essential distinction between the augmentation idea we use 
to quantize non-Lagrangian theories and somewhat similar concept of "auxiliary fields" 
jlOj . The fields are usually understood as auxiliary when they are introduced to extend 
the dynamics in such a way that the extended classical theory remains fully equivalent to 
the original one. In particular this means that the number of independent initial data for 
Cauchy problem remains the same as in the original theory. In contrast to introducing the 
auxiliary fields, the augmentation procedure results in a theory that has more degrees of 
freedom than the original one. The extra dynamics are eliminated, however, by imposing 
zero initial and/or boundary conditions on the augmentation fields. At quantum level, 
these conditions provide the absence of the "augmenting particles" in in- and out-states of 
the quantum system. It might be relevant to mention that no regular procedure is known 
yet for introducing auxiliary fields in such a way as to convert any non-Lagrangian theory 
into an equivalent Lagrangian one. In some specific models the way of introducing auxiliary 
fields is known, although it often happens that the restrictions are to be imposed strongly 
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limiting the admissible form of equations of motioij^. In contrast, the augmentation is 
a regular procedure which always works well, given equations of motion and Lagrange 
anchor. 

The paper is organized as follows. To make the paper self-contained, we review some 
recent developments in path-integral quantization of non-Lagrangian theories that includes 
basic definitions and some relevant statements from [7], [8]. In Section 2, we set up notation 
and explain some basic facts concerning the general structure of not necessarily Lagrangian 
gauge systems. We recall the notion of Lagrange structure, which contains the Lagrange 
anchor as a key ingredient, and put it in the context of Soo-algebras. The corresponding 
Subsection 2.4 is addressed to the readership familiar with basics of strongly homotopy 
algebras, others may just omit this subsection. The paper can be further understood 
without knowing the concept of S'oo-algebras, although this concept provides a natural 
homological insight into the quantization problem of non-Lagrangian dynamics. In Section 
3, we describe a BRST complex, which can be assigned to any (non-) Lagrangian gauge 
system. As input data, this complex involves the original equations of motion, generators 
of gauge identities and gauge symmetries, and the Lagrange anchor. At first, we define the 
ambient Poisson manifold that hosts this BRST complex and construct the BRST charge 
by homological perturbation theory. Further, we show that the BRST cohomology classes 
precisely correspond to the physical observables of the original (non-) Lagrangian theory. We 
also give an important interpretation of this BRST complex as that resulting from the BFV- 
BRST quantization of some constrained Hamiltonian system on the phase space of original 
fields and their sources. Section 4 is devoted to the quantization of the BRST complex. 
Quantizing the ambient Poisson manifold, we define the quantum BRST cohomology and 
present the generalized Schwinger- Dyson equation for the partition function of a (non- 
) Lagrangian gauge theory. This equation is shown to have a unique solution, which can be 
written down in a closed path-integral form. 

Section 5 contains the main results of the paper. Namely, in Section 5.1. we define 
the augmented BRST complex, which is build on the original BRST complex and carries 
all the information about the augmented theory. In Section 5.2. we unfold the structure 
of the augmented BRST charge by interpreting it in terms of equations of motion, gauge 
symmetry and Noether identity generators. As the Lagrange anchor behind the augmented 

^For example, in higher-spin field theories, the auxiliary fields can be introduced converting a non- 
Lagrangian model into Lagrangian one unless no interaction has been switched on, even though the con- 
sistent equations of motion with interaction are know for many years, see for a review and further 
references 
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theory is always nondegenerate (whatever the original anchor), the augmented partition 
function has the standard Feynman's form. Moreover, the corresponding action functional 
is proved to possess the property of space-time locality provided the original equations of 
motion and the Lagrange anchor do so. Finally, in Section 5.3. we present an alternative 
path-integral representation for the quantum averages of the original physical observables 
in terms of the augmented action functional. 

In Section 6, we apply the augmentation procedure to quantize two non-Lagrangian field 
theories: Maxwell electrodynamics with monopoles and self-dual p-forms. These models 
are known to admit no Lagrangian formulation. However, we have found non-trivial (de- 
generate) Lagrange anchors for these theories. Making use of the anchors, we apply the 
augmentation method to construct manifestly Poincare invariant partition functions for 
both the models. 



2.1. Classical dynamics. In field theory one usually deals with the space of all smooth 
maps from a space-time manifold X to a target manifold Y. The atlases of coordinate charts 
on X and Y define then a natural atlas on Y^ such that each map x : X Y is specified 
locally by a set of smooth fields x\ the coordinates on the infinite-dimensional manifold Y^ . 
Hereafter we use De Witt's condensed notation [4J, whereby the superindex "z" comprises 
both the local coordinates on X and the discrete indices labelling the components of the 
field X. As usual, the superindex repeated implies summation over the discrete indices 
and integration over the space-time coordinates w.r.t. an appropriate measure on X. The 
partial derivatives di = d/dx'^ are understood as variational ones. 

In the context of local field theory, the space Y^ is known as the space of all histories 
and the true histories are specified by a set of PDE's 



Here we do not assume the field equations to come from the least action principle, hence 
the indices i and a, labelling the fields and equations, may run through completely different 
sets. In the case where X is a manifold with boundary, Eqs. ([1]) are also supplemented with 
a suitable set of boundary conditions. Usually, the boundary conditions specify the values 
of fields and/or their derivatives up to some fixed order. Varying these values, collectively 
called the boundary data, one gets a family of different solutions to Eqs. ([1]). 

For our purposes it is convenient to think of T = {Ta{x)} as a section of some vector 
bundle S ^ M over subspace of all fields M C Y^ with given boundary data. Then the 
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Ta{x)=0. 
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set of all true histories S belonging to M is identified with zero locus of T G r(£^): 

(2) S = {x e M| T(x) = }. 

Using the physical terminology, we refer to S as the shell. Under the standard regularity 
conditions [6] , E C M is a smooth submanifold associated with an orbit of gauge symmetry 
transformations (see Eq.(l9]) below); in the absence of gauge symmetries the shell S is a 
single point of M. In the following we will always assume S to be a connected submanifold 
for each choice of boundary data. 

Thus, the classical dynamics are completely specified by a section T of some vector bundle 
£ ^ M over the space of all histories subject to boundary conditions. For this reason we 
call £ the dynamics bundle. 

2.2. Regularity conditions. To avoid pathological examples, some regularity conditions 
are usually imposed on a classical system. To formulate these conditions in an explicitly 
covariant way, let us introduce an arbitrary connection V on £" and define the section 

(3) J = WT eV{T*M ®£) 



This section, in turn, defines the M-bundle morphism 

(4) J:TM -^£, 
which is not necessarily of constant rank. 

Definition 2.1. A classical system {£,T) is said to be regular of type {m,n), if there exists 
a finite sequence of vector bundles £k ^ M and M-bundle morphisms 

(5) ^ £_^^ > £_i -^TM £ ^ £i^ ^ 

satisfying conditions: 

(a) there is a tubular neighbourhood U C M of S such that all the morphisms ([5]) have 
constant ranks over U; 

(b) upon restriction to S, the chain Q makes an exact sequence. 

This definition has several important corollaries elucidating its meaning: 
Corollary 2.1. The shell G M is a smooth submanifold with TS = Imi?|s. 



"'To simplify notation, wc will not distinguish between an A/-bundlc morphism H : £ ^ £' , the induced 
homomorphism r{H) : T{£) ^{£') on sections, and the associated section H G r{£* (g) £') ~ Mor{£, £'), 
denoting all these maps by one and the same letter H. 
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Corollary 2.2. For any vector bundle V —>■ M and a section K G r(V) vanishing on 
S C M , there is a smooth section W G T{£* (g> V) such that 

(6) K = {T,W), 

where the triangle brackets denote contraction of T and W . Informally speaking, any on- 
shell vanishing section is proportional to T . 

Corollary 2.3. When exist, the morphisms ^ are not unique off shell. Thinking of these 
morphisms as the sections of the corresponding vector bundles, one can add to them any 
sections vanishing on S, leaving the properties (a),(b) unaffected. In particular, by making 
a shift 

(7) Z^Z + Zo, ^o|s = 0, 

if necessary, we can always assume that T G kei Z . In view of the previous remark, the 
section Zq is proportional to T. 

Corollary 2.4. In the definition above we can pass from the sequence ^ to the transpose 
one by replacing each vector bundle with its dual and inverting all the arrows. The transpose 
sequence will meet the same conditions (a),(b) as the original one. 

In this paper we deal mostly with the quantization of regular (1, l)-type Lagrange struc- 
tures associated to the four-term sequences 

(8) 

The on-shell exactness at TM suggests that for any vector field V G T{TM) obeying 
condition VyTls = there exists a section e G r(jF) such that V = R{e). Combining this 
with Corollary \2.2\ we can write 

(9) K^^Ta = Uln 

for some U G T{£ ^S* ®J^). Here indices a, i, a label the components of the corresponding 
sections w.r.t. to some frames {e"} G T{J^\u), {e"} G T{£\u), and {di} G T(TU) associated 
with a trivializing coordinate chart U C M. Let {e"^} be a frame in Q over U. In view of 
Corollary 12.31 the on-shell exactness at term S implies then 

(10) Z%Ta = 0, 

if Z was chosen in an appropriate way. Relations (Q and ffTOl) have a straightforward 
interpretation in terms of constrained dynamics [6j|: the homomorphism R is identified 
with an irreducible set of gauge symmetry generators for the classical equations of motion 



10 S.L. LYAKHOVICH AND A.A. SHARAPOV 

T = 0, while the homomorphism Z generates a set of independent Noether's identities. 
Having in mind this interpretation, we term T and Q the gauge algebra bundle and the 
Noether identity bundle, respectively. Notice that the irreducibility of the gauge symmetry 
generators is provided by the on-shell exactness of (IHl) at JF, while irreducibility of Noether's 
identity generators follows from the on-shell exactness of the transpose of ^ dX Q* . 

This interpretation of homomorphisms R and Z applies to the general regular systems of 
type {m,n), except that the bases of the gauge algebra and Noether's identity generators 
may be overcomplete (reducible). A general (n + l,m + l)-type gauge theory with n > 
and/or m > corresponds to the case of ra-times reducible generators of gauge transfor- 
mations and/or m-times reducible generators of Noether's identities. The theories of type 
(0, 0) are described by linearly independent equations of motion having a unique solution. 

2.3. Lagrange structure. In the context of covariant path-integral quantization, the pas- 
sage from classical to quantum theory involves, besides classical equations of motion, one 
more geometric ingredient called the Lagrange structure [7]. 

Definition 2.2. Given a classical system {S,T), a Lagrange structure is an M-linear map 
ds : r(A"'£^) r(A"+^£^) obeying two conditions: 

(i) d£ is a derivation of degree 1, i.e., 

deiA AB) = d£AAB+ A dsB , 

for any A E r(A"f ) and B E T{A'£); 

(ii) deT = 0. 

Here we identify T{A^S) with C°°{M). 

Due to the Leibnitz rule (i), in each trivializing chart U G M the operator d^ is completely 
specified by its action on the coordinate functions and the basis sections e"" of S\u: 

(11) dsx' = V^{x)e'' , dee'' = \ci{x)e^ A . 

Applying dg to the section T = TaC"", one can see that the property (ii) is equivalent to the 
following structure relations: 

(12) dsT = {V^dm - C:,T,)e'^ A = . 

The first relation in (fTTj) means also that dg defines a bundle homomorphism V : S* —* TM. 
The section V eT{S ® TM) is called the Lagrange anchor. 
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Definition 2.3. A Lagrange structure T, ds) is said to be regular at p G M, if tliere 
exists a vicinity U (Z M of p such tliat tlie M-bundle morpliism 



(13) R®V : S^i®S* ^TM 

lias a constant ranly r over U. Tlie number r is called the rank of Lagrange structure at 
p G M. The Lagrange structure is said to be complete at p, if the homomorphism (fT3|) is 
surjective on ?7. Finally, we say that the Lagrange structure is regular (or complete), if it 
is regular (or complete) at any point of M. 



Remark 2.1. In view of Definition 12.11 the regularity of the Lagrange structure at p G S is 
equivalent to the regularity at p of the anchor morphism V : S* ^ TM, i.e., there exists a 
sufficiently small vicinity [/ C M of p G S such that V has constant rank over U. 

Remark 2.2. In the context of quasiclassical quantization we will deal with in sequel, it is 
also appropriate to introduce the notions of weakly regular and weakly complete Lagrange 
structures by requiring regularity and completeness only for the points of S. 

Theorem 2.5 ( Splitting theorem [7|). Let p E M be a regular point of the Lagrange 
structure {S,T,d£), then there is a coordinate system {y^, ...,y^ , , z^) centered at p 
together with a set of local functions S{y), E^{y), E''{y) such that equations Ta{y, z) = 
are equivalent to 

0, z' = E'iy), 



and the Lagrange anchor V = {V"^ , Vj) is given by the abelian vector distribution 

v' = o, v^, = A + ^A. 

' dy^ dy^ dz-^ 

Here the number r is the rank of the Lagrange structure at p E M. 

In case r < dimM, it is natural to call S{y) a partial action. 

Although the theorem above ensures the split of local coordinates into "Lagrangian" y's 
and "non-Lagrangian" 2;'s, it is by no means necessary to explicitly perform this splitting 
in order to develop the theory further. The subsequent formulas do not involve such a 
split. Moreover, the method is insensitive to the rank of the Lagrange anchor producing a 
well-defined path-integral quantization in the irregular case as well. 



■^Of course, in the context of infinite-dimensional manifolds the notion of rank needs clarification. An 
appropriate definition can be done, for example, in the case of local field theories. 
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Example: Let us illustrate the definitions above by an example of a Lagrangian gauge 
theory with action S{x). The equations of motion read 

(14) T = dS{x)=0, 

so that the dynamics bundle S is given by the cotangent bundle T*M of the space of all 
histories. The canonical Lagrange structure, resulting in standard quantization, is given 
by the exterior differential d : T{A''T*M) r(A"+^T*M). The defining condition for the 
Lagrange structure (fT2l) takes the form 

(15) dT = d^S = 0. 

The Lagrange anchor is defined by the identical homomorphism 

(16) V = id: TM -^TM , 

and hence the Lagrange structure is regular and complete. Suppose the action S is gauge in- 
variant. Then there exist a set of gauge algebra generators defining an M-bundle morphism 
R: ^ TM such that 

(17) {R{e),dS) = 

for any gauge parameter e G r(jF). So, equations ( fT4l) appear to be linearly dependent. 
Differentiating the last identity w.r.t. some connection V on JF (g) TM, we arrive at Rel. 
m with Ui^ = V.Rl 

Thus, we see that for ordinary Lagrangian gauge theories the dynamics bundle coincides 
with the cotangent bundle {S = T*M), the Noether identity bundle coincides with the 
gauge algebra bundle {J-' = Q), and the generators of gauge symmetry coincide with the 
generators of Noether's identities {R = Z). For a general regular system of type (1,1) 
neither of these coincidences should necessarily occur. For instance, it is possible to have 
gauge invariant, but linearly independent equations of motion; and conversely, a theory 
may have linearly dependent equations of motion without gauge symmetry. 

2.4. Soo-algebras. Recall that the conventional BV formalism for a Lagrangian gauge the- 
ory starts by introducing ghost fields to every gauge symmetry, and then an antifield for 
every field. The space of all fields and antifields is endowed with the canonical odd Poisson 
bracket ( ■ , ■ ) and the original action functional is extended to the master action S defined 
as a proper solution to the classical master equation {S,S) = 0. The classical BRST dif- 
ferential Q = {S, •), being a nilpotent derivation of the odd Poisson algebra of functions, 
incorporates then both the dynamical equations and the gauge algebra structure. Thus, the 
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odd Poisson geometry provide a natural framework for the BV field-antifield formalisir0. 
In [T], we have shown that the quantization of general non-Lagrangian gauge theories call 
for a strongly homotopical version of the odd Poisson algebras. 

Definition 2.4. An ^oo-algebra (5* for Schouten) is a Z2-graded, super commutative, and 
associative algebra A endowed with a sequence of odd linear maps Sn '■ A®"' A such that 

(a) 5„(..., ak, afc+i, ...) = a^+u •••), 
e(a) being the parity of a homogeneous element a E A. 

(b) a ^ Sn{ai, a„_i, a) is a derivation of A of the parity 
1 + Efc=J e(afc) (mod 2). 

(c) For all n > 0, 

(-l)^5'/+i(S'A,.(a<^(i), a^(fe)), 0^(^+1), a^(fc+/)) = , 

k+l=n shuflc 

where (— 1)*^ is the natural sign prescribed by the sign rule for permutation of ho- 
mogeneous elements oi, a„ G A. 

Recall that a (fc, /)-shuffle is a permutation of indices 1, 2, ...,k + I satisfying a{l) < ■ ■ ■ < 
a{k) and (7{k + 1) < ■ ■ ■ < a{k + /). 

When 5*0 = we speak about a flat ^oo-algebra. In that case 5*1 : A — > A is a differential 
with (5*1)^ = 0, and S2 induces an odd Poisson structure on the cohomology of ^i. An odd 
Poisson algebra can be regarded as an S'oo-algebra with bracket S2 '■ A ® A A and all 
other Sk = 0. In fact, properties (a) and (c) characterize Loo-algebras. We refer to [12j for 
a recent discussion of Soo-algebras. 

It turns out that any Lagrange structure of type (m, n) gives rise to a flat S'oo-algebra on 
the supercommutative algebra of sections 

m n 

(18) A = T{^'£®^ S'(n^^_fc) ® (g) S'iU'^^Ei)) . 

k=l 1=1 

Here S* stands for symmetric tensor powers (in the Z2-graded sense) and 11 denotes the 
parity reversion operation, i.e., US is a vector bundle over M whose fibers are odd linear 
spaces. By definition, 11^ = id and S'{Il£) = A*£. 



^In the physical hterature the odd Poisson manifolds are usually called anti-Poisson manifolds. Corre- 
spondingly, the odd Poisson brackets are referred to as antibrackets. On the other hand, in mathematics 
the odd Poisson algebras (brackets) are also known under the names of Schouten or Gerstenhaber algebras 
(brackets) . 
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In the next section, applying the machinery of BRST theory, we give an exphcit de- 
scription for S'oo-algebras associated with (1, l)-type Lagrange structures. Extension to the 
general Lagrange structures is straightforward. 

3. BRST COMPLEX 

3.1. An ambient symplectic supermanifold. Let {S, T, ds) be a regular Lagrange struc- 
ture corresponding to the four-term sequence ([H]). Following the general line of ideas of 
BRST theory, we realize M - the space of all histories - as the body of a graded superman- 
ifold M. The latter is chosen to be the total space of the following graded vector bundle 
over M: 

(19) Jl{T®r)®T*M®Ti{£®£*)®{g®g*). 

Here JF, and Q are the bundles of gauge algebra, dynamical equations and the Noether 
identities, respectively. The base M is imbedded into (fT9l) as the zero section. In addition 
to the Grassman parity the fibers of ( fT9i) are graded by ghost number valued in integers. To 
avoid cumbersome sign factors, we will assume the base M to be an ordinary (even) manifold 
that corresponds to the case of gauge systems without fermionic degrees of freedom. Then 
the Grassman parities of fibers correlate with their ghost numbers in a rather simple way: 
the even coordinates have even ghost numbers, while the odd coordinates have odd ghost 
numbers. The supermanifold M is also endowed with an N-grading called the momentum 
degree (or m-degree for short). 

It is convenient to arrange the information about all the aforementioned gradings of local 
coordinates in a single table: 



base and fibers 


M 


T*M 






s 


s* 


Q 




local coordinates 




Xj 






Va 


ff 


U 




e=Grassman parity 








1 


1 


1 


1 








gh = ghost number 








1 


-1 


-1 


1 


-2 


2 


Deg = momentum degree 





1 





1 





1 





1 



Table 1 



Upon splitting all the coordinates into the "position coordinates" = {x\ 0°" , rja, ^a) and 
"momenta" = (xj, c^, f/", ^^) the assignment of gradings becomes easy to see 

gh(^/) = -gHif^) , e(^/) = e{ip^) , 

(20) 

Deg(^,) = l, Deg(vpO = 0. 
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Let us denote by C°°(A/') the super commutative algebra of "smooth functions" on M . 
By definition, the generic element of C°°{N') is given by a formal power series in the fiber 
coordinates with coefficients in C°°{M). 

Fixing a linear connection V = Vjf © © Vg on © £^ © ^, we endow J\f with the 
exact symplectic structure 

(21) 00 = d{xdx' + c„Vc" + r]''Vr]a + ^^V^a) , 
where 

(22) Vc" = dc" + dx'Tff^c^ , 

and similar expressions are assumed for covariant differentials of r/'s and .^'s. Thus, C°°(A/') 
becomes a Poisson algebra; the nonvanishing Poisson brackets of local coordinates are given 
by 

{V\ Va} = , {^^, Va] = Tlf^l. , {x^, //^ = -F^,//'^ , 

{cq,, c^} = 5^ , {aSj, c"} = F^c^ , {xj, c^} = —T%Ca , 

(23) 

{e^eB} = '5^, {x„a} = rf^eB, {x.,e^} = -F,4,e^, 

Here the structure functions determining the Poisson brackets of Xi and Xj are just the 
components of the curvature tensor of V. 

Notice that the equations = define the Lagrangian submanifold 

(24) c = ii{j^ ®s)®g cU , 

and the supercommutative algebra of functions C°°(£) is naturally isomorphic to the alge- 
bra ffTSl) with m = n = 1. 



3.2. BRST charge. It turns out that all the ingredients of a classical gauge system as well 
as a Lagrange structure can be naturally incorporated into a single object fl, the classical 
BRST charge^. By definition [7j, the BRST charge is an element of the Poisson algebra 
C°^{^^) such that 

(i) e(fi) = l, gh(n) = l, Deg(f])>0; 

(ii) Q = f]^Ta + c"i?j,x, + e^Z^r/, + f^^V^x^ + ■ ■ ■ ; 

(iii) = 0. 



^Relevance of this terminology is explained in the next subsection. 
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The dots in (ii) refer to the terms which are at least hnear in rja and Cq or at least quadratic in 
Xi. Equation (iii) is known as the (classical) master equation. Conditions (i)-(iii) determine 

17 up to a canonical transformation of C°°(A/'). The existence of Q is proved by standard 
tools of homological perturbation theory [7] . 

It is instructive to consider expansion of fl in powers of momenta. In view of (i) the 
expansion starts with terms linear in (f, i.e., 



(25) n = Y,^k, Beg{nk) = k. 

k=l 

On substituting ( l25ll into the master equation (iii), we get 

(26) {ni,ni} = o, {ni,n2} = o, {^2,^^21 = -2{i^2,^^3}, etc. 

We see that the leading term Qi = VL^(pi gives rise to the homological vector field on £, 

which carries all the information about the classical system itself, with no regard to the 
Lagrange structurqj. Evaluating the nilpotency condition = to lowest order in fiber 
coordinates, one immediately recovers Rels.([9l [10]) characterizing T = as a set of gauge 
invariant and linearly dependent equations of motion, with R and Z being the generators 
of gauge transformations and Noether identities, respectively. 

The Lagrange anchor V : £* ^ TM defining the Lagrange structure for the classical 
system ( 1271) enters the next term 

(28) = = fTvix, + • • • . 

Relations fl26l) characterize ^2 as a weak anti-Poisson structure on £, i.e., Q-invariant, odd 
bivector field satisfying the Jacobi identity up to homotopy. The corresponding "weak" 
antibracket reads 

(29) (a,6) = {{fi2,a},&}, a,6GC°°(/:). 
Examining the Jacobi identity for this bracket, one finds 

(a,(6,c)) + (-l)^(^)^(^)((a,c),6) + (-1)<«)W^)+^W)((6, c), a) = 

(30) -53(Qa, 6, c) - (-l)<'^)<'')53(a, Qh, c) - (-l)(^W+^W)^W53(a, 6, Qc) 

-QS3{a,b, c) , 



'^In the usual BV theory the operator Q is known as the classical BRST differential §8.5]. 
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where we have introduced the following notation: 

(31) Sn{ai,a2,...,an) = {...{{^n,ai}a2},...,an}, akEC°°{C). 

Evidently, the weak antibracket fl29l) induces a genuine antibracket in the Q-cohomology. 

It is Rel. fl3Tl) that defines the aforementioned ^oo-structure on the supercommutative 
algebra C°°{C): By definition, each Sn is a symmetric multi-differentiation of C°°(£) and 
the generalized Jacobi identities for the collection of maps {Sn} readily follow from the 
master equation {Q, Q} = for the BRST charge. Since Deg > 0, this S'oo-algebra is 
fiat. 

3.3. Hamiltonian interpretation. In the conventional BFV approach, the BRST charge 
arises as a tool for quantizing first-class constrained Hamiltonian systems. A glance at Table 
1 is enough to see that the spectrum of ghost numbers corresponds to that of the BFV- 
BRST formalism for a first-class constrained Hamiltonian system with linearly dependent 
constraints |6]. In order to make this interpretation more explicit, let us combine the local 
coordinates with ghost numbers 1 and —1 into the ghost coordinates = (f/",c") and 
ghost momenta Vj = {ria,Co,), respectively. In this notation the above BRST charge fl25l) 
can be rewritten as 

(32) n = C'Qjix, x) + VjE'^ix, x)^ + ^PxUfjix, x)C'C' + o{V\ e) , 

where the expansion coefficients 9/ = (Ta,Ra) and = {Z'^,Z%), playing the role of 
first-class constraints and their null- vectors, are given by the formal power series in x's: 

fa{x, x) = Ta{x) + V:{x)xi + o{x^) , 

(33) Ra,{x,x) = Rl{x)xi + o{x'^) , 

Z\[x,x) = Z%{x) + o{x). 

To lowest order in C's, Eqs. ( 126|) reproduce the standard involution relations for a set of 
reducible first-class constraints w.r.t. the canonical Poisson bracket on T*M: 

(34) {ej,ej} = ufjeK, Hie, = o. 

From the regularity condition it follows immediately that the number of independent first- 
class constraints 0/ ~ is equal to dimM. In physical terms, one can interpret this fact 
concluding that the Hamiltonian system under consideration has no physical degrees of 
freedom. From the geometrical viewpoint, this implies that the equations 0/ = define a 
Lagrangian submanifold L C T*M; more accurately, L is a formal Lagrangian submanifold 
as we are not concerned with convergence of the formal series fl33|) . 
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One can also regard the constraints 0/ ~ as a formal deformation of those given 
by the leading terms of expansions ( l33l) in the "direction" of the Lagrange anchor V. 
From this standpoint, the Lagrange structure is just the infinitesimal of deformation of the 
Lagrangian submanifold Lq C T*M defined by the "bare" first-class constraints Ta{x) ~ 
and B}^{x)xi ~ 0. 

Associated with the first-class constraints 0/ ~ is the Hamiltonian action on the 
cotangent bundle of the space of all histories 

/■t2 

(35) S[\x,x\= I dt{xix' - X^Qj{x,x)) . 

The action describes a pure topological field theory having no physical evolution w.r.t. to 
t. It should be emphasized, that the "time" t is an auxiliary {d + l)-st dimension, which 
has nothing to do with the evolution parameter in the (differential) equations of motion 
Ta = 0. The true physical time is among the original d dimensions. 

The model (1351) is invariant under the standard gauge transformations generated by the 
first-class constraints and their null- vectors (1331): 



S^x' = {x\ Qi}e^ , 5eXi = {xi, <di}e^ , 

(36) 

6,X' = e' - X^'Uj.je-' + S^e^ . 

Here = (e",^") and e"^ are infinitesimal gauge parameters, and the structure functions 
U^Kji<P) are defined by (|3l. 

Imposing the zero boundary conditions on the momenta 

(37) Xi{ti) = Xi{t2) = , 

one can see [8] that the classical dynamics of the model (l35i) are equivalent to those described 
by the original (non-) Lagrangian equations = 0. 

Example: Given the Lagrangian equations of motion ( fT4l) . Lagrange anchor ( fT6l) . and 
gauge symmetry generators (|T71) . we have the following set of first-class constraints on the 
phase space of fields and sources: 

(38) fi = diS + Xi, R^ = ir^Xi . 
From the definition of gauge algebra it readily follows that 

(39) {T„T,} = 0, {R^,Rp} = U2,R, + Kpf,, {R^,f,} = Ui,f,, 
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where U^^ = diR^^ and f/^^ = XjWl^p{x). Evidently, the constraints (l38ll are reducible, 

(40) R^ = R^fi , 

and we can take {Tj} as a complete set of independent first-class constraints. The corre- 
sponding Hamiltonian action ( !35l) on the phase space of fields and sources reads 

(41) Sh[x,x,X]= [ dt{xix' - X\diS{x) +Xi)) . 

Jti 

Excluding the momenta from this action by means of equations of motion 6S/6X^ = 0, we 
obtain 

(42) Sh[x] = [ ' dtx'diS{x) = S{x{t2)) - S{x{ti)) . 

Jti 

The latter action describes two copies of the original Lagrangian theory corresponding to 
the ends of the "time" interval [^1,^2]- As there is no coupling between the fields x{ti) and 
x{t2), one can consistently restrict dynamics to either subsystem with action ±S'[a;]. This 
proves classical equivalence of the topological theory with action fHT]) to (the two copies of) 
the Lagrangian theory with action S[x]. 

3.4. Physical observables. The Poisson action of on A/" makes the space C°°{J\f) into 
a cochain complex graded by ghost number: A function A is said to be BRST-closed if 
{Q, A} = and BRST-exact ii A = {Q, B} for some B. Let H^{Q) denote the corresponding 
cohomology groups. As usual, the space of physical observables is identified with the group 
iJ°(fi), BRST cohomology at ghost number zero. 

It can be shown f7\ that the cohomology class of any BRST cocycle A with ghost number 
zero is completely determined by its restriction to M, i.e., by the function A = A\m, and a 
function O G C°°{M) is the restriction of some BRST cocycle iff 

(43) (i?(e),rfO)|s = 0, V£Gr(J^). 

The trivial BRST cocycles are precisely those for which 0|s = 0. Thus, to any on-shell 
gauge- invariant function O G C°°(M) one can associate a BRST cocycle and vice versa. 
Let [A] G H^{Q) and xq G S, then the map 

(44) [A] ^ (A) = A{xo) G R 

establishes the isomorphism H^{Q) c:^ R. Since S C M is a connected submanifold and the 
distribution R acts on E transitively (see Corollary 12. ip . the map (jHj) does not depend on 
the choice of Xq G S. By definition, (A) is the classical expectation value of the physical 
observable A. 
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4. Quantization 

In previous sections, we have described the procedure that assigns a BRST complex to 
any dynamical system, be it Lagrangian or not. The input data needed for constructing such 
a complex are the classical equations of motion and the Lagrange structure. This BRST 
complex has a clear physical interpretation as that resulting from the BFV-BRST quantiza- 
tion of the topological sigma- model (|35l) . whose target space is the cotangent bundle of the 
space of all histories. By construction, the classical dynamics of this effective topological 
theory are equivalent to the original ones for any choice of the Lagrange structure. Quantiz- 
ing now the model ( l35i) by the usual BFV-BRST method, we induce some quantization of 
the original (non-) Lagrangian theory; in so doing, different Lagrange structures may result 
in different quantizations of one and the same classical model. 

Below, we start applying the standard prescriptions of the BFV-BRST operator quantiza- 
tion to the constrained Hamiltonian system (1351) . What remains to specify is a convenient 
representation. Here we prefer to work in the coordinate (Schrodinger) representation, 
whereby a quantum state is described by a wave-function on the ghost-extended space of 
all histories. Then a physical wave-function is nothing but the probability amplitude to 
find a system developing according to a given history. For the Lagrangian systems, this 
amplitude is simply given by the exponential of the action functional multiplied by i/h. 
In the non-Lagrangian case, however, it may be a more general distribution, whose form 
strongly depends on the choice of a Lagrange anchor, (see examples in Sec. 6). 

A consistent consideration of physical states in the coordinate representation is known 
to require further enlargement of the extended phase space by the so-called nonminimal 
variables [6j. These do not actually change the physical content of the theory as one gauges 
them out by adding appropriate terms to the original BRST charge. The nonminimal 
sector just serves to bring the physical states to the ghost-number zero subspace where one 
can endow them with a well-defined inner product. We will not dwell on that in details, 
referring to the textbook [6]. From now on, VL will stand for the total (i.e., nonminimal) 
BRST charge and the phase space M will include both minimal and nonminimal variables. 

4.1. Quantum BRST cohomology. Upon canonical quantization each function on M 
turns to a linear operator acting in a complex Hilbert space Ti: 

(45) C^{J\f)3 F ^ Fe End{n) . 

A crucial step in the operator BFV-BRST quantization [6] is assigning a nilpotent oper- 
ator Cl to the classical BRST charge (l32l) . The quantum symbol of the BRST operator Cl 



21 

is supposed to have the form 

oo 

(46) Q{^,^,h) = Y,fi'^'^'\'P,'P), 

k=0 

where the leading term il^^^ is given by fl32l) and the higher orders in h are determined from 
the requirements of hermiticity and nilpotency: 

(47) Q'^ = Q, Q'^ = 0. 

It may well happen that no Cl exists satisfying these two conditions, in which case one 
speaks about quantum anomalies. In what follows we assume our theory to be anomaly 
free so that both equations (H71) hold true. 

In addition to the nilpotent BRST charge, the full BRST algebra involves also the anti- 
Hermitian ghost-number operator Q such that [Q, F] = gh.{F)F for any homogeneous F. 
In particular, 

(48) [Q,n] = Q, G^ = -g. 

Given the BRST algebra (l47l) . ( l48l) one has two BRST complexes. 

The first one is given by the space of quantum state Ti with Vt playing the role of 
coboundary operator. Under certain assumptions [6] the space Ti splits as a sum 7i = 
®n&'H'^ of eigenspaces of Q with definite real ghost number. Then : ?i" — > 7Y"^^ is the 
cochain complex of quantum states graded by ghost number. By H^^^{Q) we denote the n-th 
group of the BRST-state cohomology. 

Associated with the BRST complex of quantum states is the complex of quantum op- 
erators End{n) = ®n&End"{n). By definition, F G End"(7^), iff ad^F = [g,F] = nP. 
The corresponding coboundary operator adj^ : End"(7i) End"^^(7i) acts by the rule 
adj^F = [VL, F]. The n-th group of the BRST-operator cohomology is denoted by H^^{Q). 

The algebra of quantum physical observables and the space of quantum physical states 
are then identified with the corresponding BRST-cohomology at ghost number zero and 
the physical djTiamics are described in terms of the ifop(fi)-module H^^{Q). 

Since the BRST charge is Hermitian, the inner product on Ti induces that on the space 
of physical states Hj^p{fl). In many interesting cases, however, the induced inner product 
appears to be ill defined and needs regularization. The most popular recipe to get a regular 
inner product is to fold the BRST-closed operator e^t^'^^ ~ 1 between a pair of BRST- 
closed states 1^2) e H: 

(49) = (^i|e^[^'^]|^2), 
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K being an appropriate gauge-fixing fermion of gliost number —1. Evidently, tlie last 
expression passes to the BRST cohomology and is independent of a particular choice of K. 
(More precisely, it depends only on the homotopy class of K in the variety of all gauge-fixing 
fermions providing finiteness of ( 149|) .) Now the quantum average of a physical observable 
[6] e H^pi^) relative to a physical [|^)] G H^t{^) is given by 

(50) (o) = \ . 

4.2. Generalized Schwinger-Dyson equations. As the Hamiltonian theory we deal 
with is topological, it might be naively expected that dime H^ti^) = 1) so that the space of 
physical states is spanned by a unique (up to equivalence) BRST-closed state |$) G 7i. This 
would be quite natural because the probability amplitude must be a unique distribution on 
the space of all histories with prescribed boundary conditions. Actually, it is not always 
the case in the BRST theory: The physical dynamics may have several copies in the BRST- 
cohomology (the i7°p(f2) -module H^t{^) is generally reducible), and choosing one of them 
amounts to imposing extra conditions on the physical states [6l §14.2.6]. A guiding principle 
here is to provide a positive-definiteness of the inner product (l49l) on a superselected physical 
space. 

To be more specific, consider quantization of the gauge system ( 123|) . ( l32l) in the case 
where A/" is a superdomain endowed with canonical Poisson brackets, that is, in formulas 
( 123|) . we just take V to be a fiat connection. Furthermore, we assume that the yx^-symbol 
of the quantum BRST charge ( H6l) satisfies condition 



(51) n{^,o,h) = o. 

This property takes place for the leading (classical) term Qq = ^{f, 'f, 0) and we require 
that it holds true with account of all quantum corrections. Then, the state |$) that is 
annihilated by all the momenta, 

(52) 0im = O, 

is annihilated by the BRST charge as well. After an appropriate polarization of the non- 
minimal sector [8j, the state |$) caries zero ghost number, and hence, defines a physical 
state. In the coordinate representation, for instance, we have $(v?) = c G C At first glance 
the amplitude $, being just a constant, has nothing to do with the original dynamics, but 
that is illusion: The state \^) has an ill-defined norm in 7i, so in order to calculate the 
quantum average of a physical observable, say 1 G H^^^Q), one has use the regularized 
inner product fj49l) . but that brings an inevitable dependence of the BRST charge. In other 
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words, the information about the original gauge system enters to the state |$) imphcitly, 
through passage to the BRST cohomology. To make this dependence more exphcit one 
should consider the BRST-dual of the state |$) (see [6l §14.5.5] for general definitions). 
The dual state looks like 

(53) \<^') = \ghosts) 

and is not in general equivalent to the state |$). In the coordinate representation the first 
factor Itp), called the matter state, is described by a wave-function on M, while the second 
factor is given by a wave-function of all other coordinates. By definition, the matter state 
is annihilated by the quantum constraints: 

(54) e,|^) = o, 

where the xx-symbols of the constraint operators G/ are given by 

(55) e,(x,x,n) = ^^^^'^'^) 



ghosts=0 

It is the state that appears as physical state in Dirac's quantization method. The 
consistency of equations ( !54l) implies the following commutation relations for the quantum 
constraints: 

(56) [Qj^Qj]=Uj<jQj, 

with G's to the right of f/'s. In view of the last relation we can regard (l54l) as a non- 
abelian generalization of the Schwinger-Dyson equation to the case of non-Lagrangian gauge 
theories. The next example justifies this interpretation. 

Example: Upon canonical quantization in the coordinate representation the independent 
first-class constraints Tj in fl38|) turn to the pairwise commuting differential operators: 

(57) T, = d,S{x) - th^ , [T,, Tj] = . 

Imposing these operators on the physical wave-function ip{x), we arrive at the well-known 
Schwinger-Dyson equation in coordinate representation 

d 



(58) 



diS{x) - ih-^^. 



ij{x) =0. 



A unique (up to an overall constant) solution to this equation is given by the Feynman 
probability amplitude on M, 

(59) ij{x) = e-i^(") . 



^Hereafter the term "ghosts" refers to ah fields with nonzero ghost number. 
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One can also quantize the constraints ( l38l) in the momentum representation, which is related 
to the coordinate one by the (functional) Fourier transform. The corresponding wave- 
function 

(60) Z{x) = j Dx^(x)e-*"'^' = j Da;e-^(^(")+"'^'^ 

is nothing but the generating functional of Green's functions with x's playing the role of 
classical sources. 

In principle, one can use any copy of a single physical state in the BRST-state cohomology 
to compute the quantum average of a physical observable [O] G H^^iVL) by formula (150|) . 
It is also possible and is particularly convenient to use the asymmetric definition for the 
quantum averages: 

{^'\d\^)K 

In [8], it was shown that all such definitions give one and the same value {O). Similarly to 
the BRST-state cohomology, the BRST-operator cohomology is essentially one-dimensional 
that allows one to establish a one-to-one correspondence between the physical states and 
physical observables. Namely, given a physical observable (9, we define the physical state 

(62) \0) = d\^). 

The latter is necessarily of the form \0) = + ri|A). Using the coordinate represen- 

tation, we can rewrite fl6Tl) as 

(63) (O) = = (const) J D^OivWM . 

The last expression enables us to treat the gauge- fixed probabihty amplitude ^'^i.^) = 
{^'\(P)k as a linear functional on the space of physical observables represented by the 
physical states 0{ip) = {f\0). 

4.3. Path-integral representation. Regarding the regulator e^'^'^' in fH9l) as the evo- 
lution operator corresponding to the BRST-trivial Hamiltonian H = [Q, K], we can imme- 
diately write down the path- integral representation for the quantum average (l50l) : 

($|C)eit"'^l|$) _ 

(64) ^ I 

= (const) J V(pV^O{(p{l))exp^ dt{^np^ - {n, K}) . 
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Here the normalization constant is chosen in such a way that (1) = 1 and integration 
extends over all fields obeying 

(65) cpj{0)=ipj{l) = 0. 

These boundary conditions follow directly from definition (l52l) of the physical state |$). 
Because of (!65|) only the (^-independent part 

(66) 0{ip) = 0{ip, (^)|(^=o = 0{x) + {ghost terms) 



of the physical observable O contributes to the path integral flMj) . It is not hard to see that 
the function (166|) obeys to (and can be determined from) the following equation: 

(67) QO = {n^,O} = 0, 

Q being the classical BRST differential ( 1271) . 

Note that (Elj) is nothing but the usual Feynman's path integral for the topological 
sigma-model with action 



(68) S[^,^] 



r dt{ifj^' -{Q,K}). 

Jo 



This can be viewed as resulting from the BFV quantization of the constrained Hamiltonian 
theory (135|) . If dimX = d, where X is the initial space-time manifold, then ( l68l) defines 
a topological field theory on the {d + l)-dimensional manifold X = X x I with boundary. 
Suppose X is an orientable manifold, then so is X = X x / and each of the two orientations 
of X induces opposite orientations on the connected components of the boundary dX = 
Xo U Xi; here Xq ~ X ~ Xi and the subscripts and 1 refer to the different orientations 
of X. 

As the model ( l68l) is purely topological, there are no physical dynamics in the bulk of X. 
Put differently, all the physical degrees of freedom, if any, are supported at the boundary 
dX = Xq U Xi, where they evolve according to the classical equations of motion = 0; 
in so doing, the dynamics on Xq and Xi are completely independent of each other. Thus, 
the action (!68l) describes two copies of the same filed-theoretical model, which defer only 
by orientation of the space-time manifold X (two parallel universes). 

This classical consideration can be further promoted to the quantum-mechanical level. 
Consider the projected kernel associated with the matter state fl5^ . It can be defined by 
the path integral [6]: 

(69) ^{xi)i'{xo) = [ DipD^eT^^^'P'^^ , 
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where the sum runs over trajectories {(p^{t),ipj{t)) subject to appropriate boundary condi- 
tions at t = 0, 1. In particular, xi = x{l), Xo = x{0); the boundary conditions for the other 
variables can be found in [8]. According to our definitions, the state describes a gauge 
invariant probabUity amplitude for a field theory on AV Then i- must play the same role 
for Xi. □ Multiplying ijj by ip, we get the right probability amplitude for the field theory on 
Xq U Xi, as there is no interaction between the fields on Xq and Xi (correlations through 
the bulk of X are completely suppressed by gauge invariance). 



Example: Let us compute the quantum average flMj) for the topological model (jHj), where 
the action S is not gauge invariant. A good gauge-fixing condition in the bulk is the 
derivative gauge 

(70) F = . 

As with any abelian gauge theory, the ghost fields are decoupled from the matter ones and 
can thus be integrated out explicitly. The result is given by 

{0)=c [ DxDxDA5[x]0(x(l))e^^«[^'^'^l 



(71) =cj Dx6[x]0{x{l)) exp ^ iS[x{l)] - S[x{0)]) 

= c' j Dx(l)0(x(l))e^^["(^)l = (0|^)(V^|1), 

where 

(72) ij{x) = e^^["] , c' = c(^|l) = c / Dx(0)e-^^["(°)1 



So, up to a normalization constant, the integral ( |7T]) gives the usual quantum average of 
an observable O in the Lagrangian theory with action S. Notice that one can arrive at the 
same result by imposing the derivative gauge on the Lagrange multiplier A* = 0. 

5. Augmentation 

In previous sections, we have formulated the quantization procedure for (non-) Lagrangian 
gauge theories, which starts with the classical equations of motion and Lagrange structure 
as input data and results in the generalized Schwinger- Dyson equation for the probability 
amplitude on the space of all histories M. We have also seen that the amplitude admits a 



^In the Lagrangian field theory, for example, the probability amplitude has the form — eTi^ , where the 
action functional is given by the integral S = L of some top form L (a Lagrangian density). Changing 
an orientation of X yields S —S, hence ijj ip. 
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simple path-integral representation in terms of a Lagrangian topological field theory in the 
space-time with one more dimension. In this section, we derive an alternative path-integral 
representation for the probability amplitude of a (non-) Lagrangian theory in terms of some 
Lagrangian model on the same space-time manifold, but augmented with extra fields. The 
configuration space of the augmented field theory is taken to be the total space of the vector 
bundle E* — ^ M, the dual to the dynamics bundle E\ in so doing, the original configuration 
space M is embedded in E* as the zero section. The augmentation procedure extends the 
original (non-) Lagrangian dynamics from M to E* in such a way that the entire system 
becomes Lagrangian. We show that, at classical level, the augmented theory is equivalent 
to the original one provided that special boundary conditions are fixed for the augmentation 
fields. At quantum level, integrating the Feynman probability amplitude on E* over the 
augmentation fields yields the probability amplitude on M. 

5.1. An augmented BRST complex. Augmentation of the original dynamics on M 
implies a consistent extension to E* of the original equations of motion, gauge symmetries, 
Noether identities, and the Lagrange structure. As a practical matter, it is convenient to 
make these extensions not at the level of the space of all histories, but augmenting the 
ambient symplectic manifold A/", which already involves all necessary ghost fields of the 
original theory. The overall result of these extensions turns out to be just a "duplication" 
of the ambient manifold. More precisely, the manifold A/", considered as the total space of 
the vector bundle f|T9|) . is replaced with A/lug = A/" © 11 (A/" ©TM). The table below contains 
the data on various gradings assigned to the fiber coordinates of Yl{M © TM): 



fibers 


E* 


E 


g* 


g 


T*M 


TM 






fiber coordinates 


ya 


Vb 




cb 


Vi 


fjj 






e=Grassman parity 








1 


1 


1 


1 








gh = ghost number 








1 


-1 


-1 


1 


-2 


2 


Deg = momentum degree 





1 





1 





1 





1 


deg = resolution degree 


1 





1 





1 





1 






Table 2 

In order to compare the ghost numbers of the new and old fields it is convenient to 
assemble the augmentation fields into "position coordinates" and "momenta" : 



(73) 



= iV^, a, Z/^ C^) , ^aug = iv\ r, Va, Ca) 
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Then we have 
(74) 

Deg((^Lj = 1 , Deg(^D = , 

and 

(75) ghi^T') = sH^i) - 1 > sH^iJ = sW) + 1 • 

So, the "duphcation" of the ambient manifold Af is accompanied with reversion of parities 
and shift of ghost numbers. 

As a next step, we extend the exact symplectic structure (pT]) on Af to that on A/lug by 
setting 

(76) ^aug =UJ + c?(<^LgV(^f ^) , 

V being some connection on A/" © TM. 

Finally, the original BRST charge f2 on A/" is extended to A/lug as 



n , 

n=Q 



(77) (^aug = n + J2^n, deg(n„) 
Here 

(78) fio = 
and the higher orders in the resolution degree are determined from the master equation 

(79) {^aug) ^aug} . 

Let us show that the last equation has a solution indeed. To this end, we introduce the 
following pair of nilpotent operators: 



d_ 

o au 



S = ViTr^^ (5^ = 0, deg((5) = -1 



(80) 

5* = ¥.rA, {S*y = 0, deg(5*) = l. 
It is straightforward to check that 

(81) 66* + 6*6 = N, N = Ni + N2, 

where the operator 
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counts the momentum degree of the "old" variables (see Table 1), while 
(83) = ^T'gir, 

counts the resolution degree of augmentation fields. If we regard 6 as the differential of the 
cochain complex C°° (A/lug), then 6* becomes a homotopy for with respect to 6. As a 
result, all the nontrivial (5-cocycles are nested in the subspace ker C C°°(A^aug)- 

Now applying the standard technique of homological perturbation theory [6], we can 
prove the following statement. 



Proposition 5.1. There is a unique BRST charge ( f?^ satisfying the master equation ( [yg| j 
and the condition 

(84) 5*(^]aug-^^-^^o) =0. 

Proof: Expanding the master equation ( 1791) with respect to the resolution degree, we arrive 
at the following sequence of equations: 

(85) 5n„+i =5„(r]o,-,^^n), neN, 
where 

n 

(86) Br. = P„ {{n, Ur,} + Y,{^n-s, ^s}) , 

s=0 

and Pn is the projector on the subspace of functions of resolution degree n. We can solve 
these equations in series starting with 6Qi = Bq. Since degi?„ = n and the function 
Bq = {Q, Qq} contains no terms of zero momentum degree w.r.t. the old variables, the 
operator is invertible on the subspace W C (A/lug) spanned by all B's. The condition 
6Bn = is then necessary and sufficient for the n-th equation (1851) to be solvable. The 
closedness of i?„ is established by induction on n, just putting successive restrictions on the 
resolution degree of the Jacobi identity {Q, {Q, fi}} = 0. 

Finally, applying the operator 6* to both sides of Eq. fl85l) and using Rels. flSTl) . fl8^ . we 
get the following recurrent relations for the homogeneous components of Q: 

(87) = 6*{N\w)-'B^{Qo, ^n) • 

Here we have used the fact that the operator 6* commutes with and, as a consequence, 
with {N\y)~^- By construction, 6*Qn = 0, Wn > 0, so that the augmented BRST charge 
flaug nieets equation (!84l) . □ 

In sequel we will need the following property of the augmented BRST charge. 
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Proposition 5.2. // the augmented BEST charge ( fTTp satisfies 1^8^, then 
(88) (^]aug-^^-^^o)|^L,=o = 0. 

Proof: by induction on resolution degree. 

5.2. Interpretation. Given the augmented BRST charge fl77|) . one may ask what is a 
classical theory this BRST charge corresponds to (or can be derived from). According to 
the general definitions of Sec. 3, the equations of motion, gauge symmetry and Noether 
identity generators, as well as the Lagrange anchor, enter the BRST charge as coefficients 
at lower powers of fiber coordinates with certain ghost numbers and momentum degrees 
(see relation (ii) at the beginning of Sec. 3. 2). So, to identify all the key ingredients of 
the underlying gauge dynamics we are going just to evaluate the appropriate terms in the 
augmented BRST charge. Thus, the equations of motion define the terms that are linear 
in fiber coordinates: 



^9) 



aug 



Taix) 



0. 



on 



aug 



V,T,y'^ + o(i/2) = 0. 



As is seen, the first group of equations coincides with the original equations of motion on 
M. The absence of ^-contributions to these equations is guaranteed by Proposition 15. 2[ 
So, the original dynamics on M are completely decoupled from the augmented system (!89|) . 
The second group of equations, being at least linear in y's, admit a trivial solution = 0, 
which can be singled out by imposing zero boundary conditions on y's. 

In general, the augmented equations of motion (l89l) are both gauge invariant and linearly 
dependent. It follows from definiens of Sec. 3. 2 that the gauge algebra generators are given 
by 



(90) 



so that 



(91) 



aug 



dc"dxi 



d d^n 



aug 



dc^dya 



dx'' 
d 



+ 



aug 



+ 



dc^dy, 



aug 



dc^dxj 



d 

^ = Zl{x] 



7, dx 



R^Tq 



d 

Qya 



+ o{y) 



for some structure functions U. The Noether identities have the following form in the 
augmented theory: 



(92) 



0, 
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where 



(93) 



aug 



"aug 



aug 



aug 



o{y) 



o{y)- 



As is seen from Rels. ( l9Up . there are two types of gauge symmetry transformations in 
the augmented theory. The first ones, generated by Rq,, are just extensions to £* of the 
original gauge symmetries. The second type transformations, generated by R^, start from 
the vertical vector fields on £* associated with the Noether identity generators Z^. 

Looking at the generators of Noether's identities ( l93l) . one can observe a mirror inversion 
in the structure of the gauge symmetry generators: The generators Ra on M give rise 
to the Noether identity generators Z^,, while the generator Z^ is just a continuation of 
corresponding Noether identity generator from the original theory. 

We thus conclude that the numbers of Noether's identities and gauge symmetries coin- 
cide in the augmented theory. Furthermore, the expansion in the augmentation fields 
starts with the same terms for both sets of generators. And one can further deduce that 
the generators of identities (p3|l coincide with the generators of gauge symmetries (!90|) . 
Such a pairing between Noether identities and gauge symmetries is characteristic for the 
Lagrangian dynamics. 

To further elucidate the meaning of the augmented BRST charge in terms of the phase 
space of fields x\ y"" and their sources Xj, we introduce the following collective notation: 



(94) 



[Xi,ya) 



Then the deformed phase-space constraints associated with the augmented equations of 
motion ( l89l) are given by 

~ on. 



"aug 



(95) 



T.(0) + 



c 



S*®T''M®£ 



■bk. 



0. 



k=2 



According to our definitions, the coefficients V^{(j)) in (1951) are to be identified with the 
components of the Lagrange anchor. Using the recurrent relations fISTI) . we find 



(96) 



+ o{y) . 
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As is seen the augmented Lagrange anchor is always nondegenerate and its inverse has the 
form 

(97) K = V-'=\ + o{y) . 

To make contact with the definitions of Sec. 2, we identify the total space of the tangent 
bundle T£* with the total space of £* © TM © £* and the total space of T*£* with that of 
£(BT*M(B£* by making use the linear connection V on ^ M. Upon these identifications, 
the bundle map T£* £* goes into the bundle map £* Q)TM Q)£* £* (projection on the 
third factor) and the same is true for the cotangent bundle T*£*. Now we can summarize 
the discussion above as follows. 

Proposition 5.3. The augmented BRST complex describes a complete Lagrange structure 
of type (1,1) associated to the on-shell exact sequence 

(98) ^ V ^ r£:* ^ T*£* ^ V ^ , 

where the gauge algebra (= Noether identity) bundle V is the vector bundle with the base 
£* , total space ® £* ® G, o,nd the bundle map p : T ® £* ® Q £* (projection on the 
second factor). 

The completeness of the augmented Lagrange structure has two immediate consequences. 
First of all, the constraints ( l95i) define a Lagrangian submanifold in the augmented phase 
space £* © T*M © so that the rest of the constraints, namely, the constraints 



(99) 



= Raix,x) + o{y) , 

€*®T*M®£ 

= Zl{x,x)ya + o{y) 

e*®T*M®e 



associated with the gauge symmetry generators fl90l) . fl33|) . are given by linear combina- 
tions of (1951) . The second consequence is that, according to Theorem 12. 5[ the augmented 
equations of motion (|89|1 are equivalent to Lagrangian ones. 

To get an explicit expression for corresponding action functional, one has just to resolve 
the first-class constraints fl95|) with respect to momenta 0^. This can always be done at 
least perturbatively. As a starting point, we rewrite the constraint equations (p5l) in the 
following equivalent form: 



(100) 0, = -a|t5-a|X;C- 



'Pai ' ' ' ) 



k=2 
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where A is defined by relation (1971) . Then, taking cj) = —AT as zero order approximation 
and iterating these equations ones and again, we finally arrive at the equivalent set of 
first-class constraints 

(101) t; = 0,-t;(0)^o, 

where 

CO 

(102) n = F'--'- (</>)T5, (0) ■ ■ ■ T-,^ {<P) = mlTM . 

k=l 

The constraints fllOip . being resolved w.r.t. momenta (pa, are to be necessarily commuting, 

(103) {n,T^J = 0, 

that amounts to existence of an action functional S{(f)) such that 

(104) T-M = d,S{<P) . 

Thus, the augmented equations of motion are equivalent to the Lagrangian equations (11041) 
with A playing the role of integrating multiplier. Finally, using the standard homotopy 
operator for the exterior differential, we can reconstruct the action as 

(105) S{(j)) = 0" / T'^is(j))ds + (const) . 



Up to the second order in y's and an inessential additive constant the action reads 

(106) S{x, y) = T,(x)y" + GMy'^y' + o{y') . 
Here the symmetric matrix 

(107) Gab = Vy^Tb + ViV^Ta 

can be thought of as a generalization of Van Vleck's matrix. It is the matrix that defines 
the form of the first quantum correction to the classical average of physical observables |7] . 
More explicitly, the equations of motion ( I104p following from variation of (1105^ read 

(108) = — = Ta{x) + o{y) = , = — = S.T.y'^ + o{y') = . 

As is seen, the dynamics on M do not decouple from those on the augmented configuration 
space S* for arbitrary boundary conditions of y's, as opposite to (1891) . Nonetheless, imposing 
zero boundary conditions on y's, we can satisfy the second group of equations in (11081) with 
y = and x is arbitrary. Then the first group of equations reduces to the original equations 
of motion on M. 
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An important observation on the action (11061) is that it has the form of local functional 
whenever the augmented constraints (l95i) are locaS. Indeed, the only place where non- 
locality could emerge is the inversion of the augmented anchor ( l96l) . But, as is seen from 
( 197|) . the inversion procedure, being performed perturbatively in y's, does not spoil locality. 
Therefore, the equivalent Lagrangian equations fll04p are local and so is the action functional 

(unsD. 

5.3. Quantizing non-Lagrangian dynamics via augmentation. As the augmented 
theory is always Lagrangian, its probability amplitude has the standard form 

oo 

(109) VI/ (x, y) = ei^'^^'y^ , Sn{x, y) = Y1 ^ y) . 

n=0 

Here the leading term So{x,y) is given by the classical action (11061) and the other terms 
can be regarded as quantum corrections to the naive path-integral measure dxdy on S*. By 
definition (Ell), the probability amplitude (I109p is a unique solution to the Schwinger-Dyson 
equations 

(110) fj^{x,y) = 

associated with the (over) complete set of augmented constraints ( l95l) . Notice that the 00- 
symbols of the quantum constraint operators in (IllOp may defer from fl95|) by some quantum 
corrections in h. These corrections can be systematically derived by solving the quantum 
master equation ^2^^g = for the augmented BRST operator. 

What we are going to show in this section is that integrating the amplitude (11091) of 
y's, we get the solution to the original Schwinger-Dyson equations fl5^ . In other words, 
averaging the augmented probability amplitude (11091) over the fibers of the vector bundle 
S* —>■ M yields the probability amplitude for the original (non-)Lagrangian dynamics on 
M. We prove this statement under the following technical assumptions: 

a) The normal symbol of the augmented BRST operator fiaug, which may differ from 
the classical BRST charge (1771) by some quantum corrections, still obeys Rel. (l88l) . 

b) Both the augmented and original constraint operators are hermitian (w.r.t. the stan- 
dard inner product associated with the translation-invariant integration measures 
dxdy and dx on £* and M, respectively). 

Note that the second condition follows from hermiticity requirement for the BRST operator 
provided that the fields f]^ are chosen to be real, i.e., [f]^)* = f}^ . 

^''i.e., given by ordinary functions of fields (0^, (^j) and their derivatives up to some finite order. 
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Proposition 5.4. Under the assumptions above, 

(i) the physical observables of the original theory are also observables of the augmented 
theory; 

(ii) the functional 

(111) i^{x) = j dy^{x,y), 

where \l/(a;, y) is the Feynman amplitude U09\) and the integral is taken over all y 's satisfying 
zero boundary conditions, obeys the original Schwinger-Dyson equations [54\ ) in coordinate 
representation; 

(Hi) let O be the physical observable associated with an on-shell gauge invariant function 
O E C°°{M) of the original theory, then the quantum average (61\) is given by 

(112) (O) = {const) j dydxO{x)^{x,y) . 

Remark 5.1. The integrals (11 11 1) and (11121) . as they stand, are well defined only for theories 
of type (0,0). In presence of gauge symmetries and/or Noether identities one should treat 
the action Sh{x, y) within the usual BV quantization method. This implies extension of the 
augmented configuration space £* hy ghost fields and imposing gauge fixing conditions that 
effectively reduces integration (11121) to the space of gauge orbits. It is a perfectly standard 
technology and we will not dwell on it here. 

Proof. Statement {Hi) is an immediate consequence of {i) and {ii). 

We start with proving {ii). By definition, the amplitude \E'(x, ?/) is the matter state 
annihilated by the operators of augmented constraints. In particular, it is annihilate by the 
constraint operators that are extensions to £* of the original constraints (133|) . We have 

(113) e,^(a;,2/) = 0, 

where the constraints 0/ = (Ta,RQ,) are defined by Rels. fl95l) and (p9l) . The condition 
(188|) . being imposed on the normal symbol of the augmented BRST charge, suggests the 
following structure for the 00-symbols of the constraint operators: 

(114) @i = Qi{x,x) + o{y). 

That is the difference 0/ — 0/ between the original and augmented constraints is not only 
at least first order in y's it is also at least first order in y's. Now multiplying (11131) on an 
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arbitrary function of and integrating the result over we get 

= / dxdy^{x)@i'^{x,y)= / dxdy {@\^{x))'${x,y) 

(115) 

= J dxQ\^{x) J dy'^{x,y) = J dx^{x)Qi J dy'^{x,y). 
Here we have used the hermiticity requirements 

(116) el = e/, el = e/, 

and ReL (11141) . Since the function $(x) is arbitrary, the identity (IllSp is equivalent to the 
desired one 

(117) e, / dy^>{x,y) = Q. 



Now, let us prove {i). As we have already mentioned in Sec J3.4l the space of physical 
observables is canonically isomorphic to the space of on-shell invariant functions on con- 
figuration space modulo trivial ones. In particular, a function O G C°°{£*) gives rise to a 
BRST invariant function on A/lug, i-e., an observable of the augmented theory, iff 

(118) R„0 = WST, + W;T,, R^O = W^T, + W^T, 

for some W's. Here the generators of the augmented gauge algebra may differ from ( l99l) by 
quantum corrections. Due to Proposition 15.21 and our assumptions these generators have 
the following structure: 

(119) R„ = i?^(x)— + <(x,y)— , R^ = R'X{x,y)- 



dx' '^'9?/'*' ''''dy 

It remains to observe that for a y-independent function 0{x), Eqs. flllSp reduce to the 
on-shell invariance condition 



(120) RaO = W:Ta. 

This completes the proof. □ 

Example: Consider a Lagrangian theory with action S{x). In this case, the dynamics 
bundle coincides with the cotangent bundle T*M of the space of all histories. For simplicity 
sake assume that T*M admits a flat connection. Given the canonical anchor (|T6l) . the 
augmented constraints ([95]) on TM © T*M © T*M read 

(121) T, = d,S{x) +x,-y,, T[ = d,S{x + y) - d,S{x) - x, . 
The action of the augmented theory (11061) takes the form 

(122) S{x, y) = S{x + y)- S{x) = y'd,S{x) + \jYd,d^S{x) + ■■■ . 
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After normalization, the quantum average of a physical observable 0{x) coincides with its 
usual value 



Of course, in the presence of gauge symmetries both these integrals are to be understood 
as integrals over the space of gauge orbits rather than over TM or M. 

In this Lagrangian case, it is quite natural to interpret the augmentation fields as the 
variations of the original fields and this interpretation is automatically consistent with 

the zero boundary conditions for |/'s. 



In this section, we demonstrate by examples what the Lagrange anchor can look like in 
non-Lagrangian relativistic field theory and how the general formalism described in previous 
sections works in practice. As the examples we consider two illustrative non-Lagrangian 
models: Maxwell electrodynamics with monopoles and self-dual p-form fields. 

The Maxwell equations are considered in terms of the strength tensor and, in this formula- 
tion, they are not Lagrangian even without magnetic currents. If a magnetic monopole was 
point-like and satisfied the Dirac quantization condition, the theory would admit an equiv- 
alent Lagrangian formulation in terms of vector potential. We consider generic magnetic 
and electric sources, so the theory does not have any Lagrangian reformulation, although 
it still has a nontrivial Lagrange anchor, that is sufficient for a consistent quantization of 
the model. The observation about the structure of the Lagrange anchor might also be 
instructive for other non-Lagrangian field theories formulated in terms of strength tensors. 

Studying the second example, we reverse the order of exposing the quantization procedure 
as compared to that described in Sections 4 and 5. Given the equations of motion for a 
non-Lagrangian field theory in d dimension and a compatible Lagrange anchor, the general 
method allows one to cquivalcntly reformulate this theory as a topological Lagrangian field 
theory in ri + 1 dimensions, with the original dynamics being localized at the boundary 
of this {d + l)-dimensional space-time. This also allows for a reverse consideration: one 
can start with an appropriate {d+ l)-dimensional, topological Lagrangian theory and then 
identify "original" d-dimensional field equations and a Lagrange anchor in the action of 
the topological theory. In practice, this can be an instructive scheme for identifying those 
non-Lagrangian models that admit Poincare covariant Lagrange anchors. To exemplify this 
idea, we take the Chern-Simons theory in 4n -|- 3 dimensions and reinterpret it as resulting 
from some quantization of self-dual {2n + l)-form fields in 4n -|- 2 dimensions. 



(123) 




6. Examples of quantizing non-Lagrangian field theories 
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6.1. Maxwell electrodynamics with monopoles. Consider the Maxwell equations with 
electric and magnetic currents: 

(124) d)F = I, d)F = J. 

Here F is the strength of the electromagnetic field, whose Hodge dual is denoted hy F = *F, 
J and / are the electric and magnetic currents, respectively, and S = *d* is the adjoint 
exterior differential. As a consequence of Eqs. (11241) . the currents J and / are conserved, 

(125) rff/ = 0, SJ = Q. 

Clearly, Equations (11241) . as they stand, are not Lagrangian, even if we set J = 0. (The 
number of equations is less than the number of fields). 

Let us introduce the source P which is canonically conjugate to the field F . With the 
field F, being a 2-form on the space-time manifold, the source P is a bivector field on the 
same manifold. The canonical symplectic structure on the cotangent bundle of the space 
of all histories is given by 

(126) u = J 5FA5P', 

where the 2-form P' is obtained from P by lowering the upper indices with the space-time 
metric. 

Consider now the following set of first-class constraints on the phase space of fields and 
sources: 

T^ = SF-I^0, T'^ = d^{F + P') - J ^0, 

(127) 

{T^T''} = 0, a, 6 =1,2. 

These constraints are obtained from (11241) by adding the momentum depending term d^P' 
to the second group of equations. It is the term that defines the canonical Lagrange anchor 
for the Maxwell electrodynamics [7j. Observe that the anchor is regular but not com- 
plete (see Definition 2.3). The physical meaning of this incompleteness can be understood 
in the following way. Let / = 0, then the first group of equations (11241) expresses the 
closedness condition for the strength form F. The absence of momentum contributions to 
the corresponding constraints ^ implies that we consider these equations as being 
pure non-Lagrangian in the sense of Theorem 2.5. Hence, no quantum fiuctuations violate 
the closedness condition dF = 0, that guarantees the existence of a local gauge potential 
A = d~^F both at classical and quantum levels. 
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Notice that the constraints (11271) are hnearly dependent, d'^T"- = 0, while the classical 
equations of motion (I124p are not gauge invariant. Thus, according to Definition 2.1, we 
have a theory of type (0, 1). 

Upon canonical quantization, the constraints (1127p turn into the following Schwinger- 
Dyson operators: 

(128) f^ = d^F-I, = d^ [F -th^) - J, 



6F\ 

F' being the contravariant strength tensor of electromagnetic field. The corresponding 
Schwinger-Dyson equation for the probability amplitude 

(129) r"^[F] = 
is satisfied by 

(130) ^[F] = A[T^]e^^[^] . 
Here 

(131) S[F] = j ^GdF A *dF -GFAdJ, A[T^] = j DC5[T^ - dC] , 

C is an auxiliary 0-form, and G is the inverse of the Laplace operator □ = dd'^+d'^d. One can 
easily see that the distribution A[T^], considered as the functional of F, is supported at the 
points where T^[F] = so that T^A[T^] = 0. (A naive solution to the last equation, namely 
A[T^] = 5[T^], is ill defined because of linear dependence of the constraints T^.) Notice 
that the amplitude (I130p is non-Feynman: it is a nearly everywhere vanishing distribution 
on the configuration space of fields rather than a smooth, complex-valued function with 
absolute value 1. This fact is a direct consequence of incompleteness of the Lagrange 
anchor discussed above. 

Passing to the momentum representation, we get the generating functional of Green's 
functions 



(132) 



Z[P] = J DFvl/[F]e^^^'^^ = e^^t^l , 
W[P] = j ^Gd^P' A*SP' - F A*P' , 



where 

(133) F = dGJ + *dGI 

is the mean electromagnetic field produced by the sources / and J. As with any free 
theory, the mean field F satisfies the classical equations of motion (11240 . One can also see 
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that the propagator {F{x)F{x')) for the field F coincides with corresponding expression 
{dA{x)dA{x')) in the Maxwell electrodynamics with action S[A] = | / dA A *dA. 

The probability amplitude fll30p can also be arrived at by applying the augmentation 
method. By definition, the augmentation fields are the sections of the bundle which is 
dual to the dynamics bundle of the theory. So we introduce the 1-forms B^, a = 1,2; the 
pairing between the equations of motion and augmentation fields is given by the integral 
J A *Ba. Since the constraints fll27p are linear in fields and momenta, the action S[F, B] 
of the augmented theory is at most quadratic in F and 5". Specializing the general formulas 
f ll06p and fll07p to the case at hand, we find 

(134) S[F, B]= j {d^F -I) A *Bi + {d^F - J) A *B2 + ^dB2 A *dB2. 

The Noether identities between the original equations of motion ( ]124p give rise to the gauge 
invariance of the action (11340 : 

(135) 6eBa = dea, a = 1,2. 

We can fix this arbitrariness by imposing the Lorentz gauges d^B'' = on the augmentation 
fields and adding these constraints to the action (I134p with the Lagrange multipliers C". 
Then the gauge-fixed action reads 



(136) 



5gf [F, B, C] = S[F,B] + J SBa A 



According to Proposition 5.4, the (non-Feynman) probability amplitude (I130p admits the 
following path-integral representation in terms of the local action fll36p : 

(137) ^[F] = {const) j DBDCe^^>^'^^'^'^^ 

Of course, in the case under consideration, one can verify the last equality directly, either by 
calculating the Gauss integrals over _B's and C's or substituting (11370 into the Schwinger- 
Dyson equation (11290 and differentiating under the integral sign. 

Given the probability amplitude (11300 . the quantum average of a physical observable O 
is defined by the path integral 

(138) (O) = j DFO[F]^[F] . 

In case / = 0, one can solve the constraint = d'^F in terms of the gauge potential A 
obeying the Lorentz gauge-fixing condition, 

(139) F = dA, d^A = 0, 
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and integrate the pre-exponential A-functional in f ll30p as 



(140) 



DCDFA[d)F + dC] 



DA5[SA] . 



Then the integral fll38p takes the form 



(141) 



(0)= / DAO\F{AmA], 



where 



(142) 




is nothing but the usual probability amplitude for the electromagnetic field subject to the 
Lorentz gauge. 

6.2. Self-dual p-form fields. It has long been known that the quantization of chiral 
bosons in (4n + 2)-dimensional space-time is closely related with the quantization of Chern- 
Simons theory in the space-time with one more dimension |14|. Roughly speaking, a phys- 
ical wave-function of Chern-Simons fields on a {An + 3)-dimensional manifold M. can be 
treated as a probability amplitude (or partition function) for the self-dual fields living on 
the boundary of M.. For a recent discussion of the relationship between self-dual fields and 
the Chern-Simons theory we refer the reader to [15]. Below we justify and reinterpret this 
ad hoc quantization technique within the general method of Sections 4 and 5. 

Our starting point is the Chern-Simons action for the {2n + l)-form field F G A^''+^{M) 
on a (4ra + 3)-dimensional manifold M. 



D being the exterior differential on M. Assume that Ai = M x I, where / = [0, 1] C M and 
M is a compact {4n + 2)-dimensional manifold without boundary. Then dAi = M U M. 

Using the product structure of the manifold one can globally decompose the field F 
and the operator D as 



Here H G A^"+^(M) and B G A^"(M) are the one-parameter families of differential forms 
labelled by t G [0, 1] and d is the exterior differential on M. In this notation, the action 
f ll43p takes a simple Hamiltonian form ( l35i) . if one identifies t with evolution parameter: 



(143) 




(144) 



F = H + B Adt, 



D = dtAdt + d, 



(145) 
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The fist term in (11451) defines (and is defined by) a symplectic structure on A^""'"^(M); the 
corresponding symplectic 2-form reads 

(146) uj= [ 5HA5H. 

The field B plays the role of the Lagrange multiplier to the first-class constraints 

(147) T = dH^O. 

It is convenient to treat T as a linear functional (de Rham's flux) on the space of 2n-forms: 

(148) = / aAdH, VaeA2"(M). 

Then one can easily check that the constraints (11471) have vanishing Poisson brackets, 

(149) {T[a],T[p]}= [ da A dp = 0, Va, /3 G A2'^(M). 

Since d"^ = 0, these constraints are reducible, dT = 0, and one can further deduce that 
the order of reducibility is 2n. 

The Hamiltonian reduction by the first-class constraints (11471) leads to a finite dimensional 
phase space. We have a rather explicit description of the reduced phase space due to the 
Hodge decomposition 

(150) A2"+1(M) = rfA2"(M) © dU^"+^(M) © A^"+^(M) . 

Here d^ : A™'(M) A'^~^(M) is the adjoint differential constructed by some Riemannian 
metric on M, and A^^^(M) is the subspace of harmonic forms. According to the Hodge 
theory, the space A^'^^(M) is naturally isomorphic to the de Rham cohomology group 
H'^"'~^^{M). The first-class constraints (11471) single out the coisotropic subspace of (i-closed 
forms, whose complementary isotropic subspace is given by the d'^'-exact forms. Notice that 
both these subspaces are Lagrangian iff H'^^~^^{M) = 0. 

The Hamiltonian flux generated by the first-class constraints changes any {2n + l)-form 
H by an exact one: 

(151) 6eH = {H,T[e]} = de. 

Taking the quotient of rf-closed (2?7,-|-l)-forms by rf-exact ones, we obtain the physical phase 
space of the model, which is apparently isomorphic to the (finite dimensional) subspace of 
harmonic forms on M: 



(152) 



A2'^+1(M)//T ~ A^"+i(M) ^ H^''+\M) 
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We are lead to conclude that the model under consideration is not topological unless 
iJ^"+^(M) = 0. (We define a topological theory as a theory without physical degrees 
of freedom.) To get rid of the physical modes and obtain a pure topological model we 
can restrict the dynamics on the affine subspace A^""*"^ (M) C A^""*"^ (M) constituted by the 
forms H = Ha + Ho, where Hq e dK^''^{M) © d)h?''^'^{M) and H^ is a time-independent 
harmonic (2n + l)-form representing the de Rham class a = [Ha] G iJ^"^^(M). This re- 
striction is compatible with dynamics. Indeed, the equations of motion following from the 
Hamiltonian action f ll45p read 

(153) H = dB, dH = 0. 

So, the de Rham class [H] of the closed form H does not change with time and can thus 
be regarded as a (topological) integral of motion. Moreover, the embedding A^"+^(M) C 
A^""'"^(M) is symplectic, i.e., the restriction of the 2- form (11461) to A'^~^^{M) is nondegen- 
erate. 

To further proceed with the interpretation and quantization of the Chern-Simons theory 
on the product manifold = M x J, let us endow M with a Lorentzian metric. (A 
necessary and sufficient condition for such a metric to exist is that the Euler characteristic 
x(M) be zero.) Then, the corresponding Hodge operator * squares to +1 on the middle 
forms so that any (2n + l)-form H admits a unique decomposition in the sum of its self-dual 
and anti-self-dual parts: 

(154) H = H+ + H- , *H^ = ±H^ . 

Since uj{6H^, 6H^) = 0, we have a natural polarization of the phase space A^"+^(M) given 
by the two complementary Lagrangian subspaces of self- and anti-self-dual forms: 

(155) A2"+1(M) = Al^+\M) © A?"+i(M) . 

Let us regard the fields H~ as the "momentum coordinates" canonically conjugate to the 
"position coordinates" H~^ and rewrite the Hamiltonian action f ll45p as 

(156) S= [ dt [ H' AH+ -BA d{H^ + H') 

J I J M 

Upon restriction to A^"+^(M) this action describes a topological field theory and its form is 
identical to the form of the topological action fl35|) . Recall that the latter was constructed 
on general grounds starting from some classical (not necessarily Lagrangian) equations of 
motion supplemented with an appropriate Lagrange structure. Proceeding now backward, 
we can readily reinterpret the topological model (11561) in terms of non-Lagrangian dynamics 
on M. 
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Namely, the classical equations of motion are to be identified with the momentum inde- 
pendent terms in the Hamiltonian constraints 



(157) T = d{H+ + H-) ^0. 

With our choice of the phase-space polarization this yields the closedness condition for the 
self-dual form H~^, 

(158) dH^ = . 

Similar to the Hamiltonian constraints (11571) . these equations are 2n-times reducible, and 
hence they define a regular gauge theory of type (0, 2n) (although there is no gauge in- 
variance in the usual sense). A non-Lagrangian nature of equations fll58p was discussed at 
length in |13j . 

The second term in the Hamiltonian constraints (11571) . namely dH~ , is linear in momenta 
and should be identified with the Lagrange anchor. Zero boundary conditions (!37l) on the 
momenta 

(159) H-\t=o^, = 

ensure the equivalence of the classical dynamics (11531) and (11581) . where both H and 
belong to A^""'"^(M). Fixing the de Rham class [H^] = a of a solution to equation (I158p is 
quite similar to fixing the boundary conditions for a field theory on a bounded space-time 
domain. 

Quantizing now the Hamiltonian constraints (I157p in the coordinate representation, we 
get the following Schwinger-Dyson operator: 

(160) f = d(^H^- ^hJ^^ . 

Here H'^ is the bivector on M obtained from by rising indices with the help of the 
Lorentzian metric. The probability amplitude on the configuration space A'^~^^{M) fl 
A^~^^{M) is determined by the equation 

(161) f^[i/+] = 0. 

We use the augmentation method to write down an explicit path-integral representation 
for \E'[if"'"]. To this end, we introduce the augmentation fields C G A^"(M) whose configu- 
ration space is dual to the linear space A^"+^(M) of the field equations (11581) . and apply the 
general formulas (11060 . (11071) . and (I109P to construct the Feynman probability amplitude 
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of the augmented theory. The result is almost obvious: 
(162) 

S[C, H+] = / --dC A *dC + H+ AdC. 
The corresponding equations of motion read 

(163) dH+ = 0, {dC)- = (}. 

The augmented theory is seen to describe the pair of self-dual fields: one in terms of the 
"strength tensor" and another one in terms of the "gauge potential" C. Notice that 
the Noether identities for the non-Lagrangian equations of motion fll58p . i.e., d{dH~^) = 0, 
reincarnate as the gauge transformations of the augmentation fields: 

(164) C^C' = C + dA, VAgA^"^^(M). 

Integrating formally the amplitude (11621) over the fields C, we obtain the probability 
amplitude for the self-dual field H~^, 

(165) ^>[H+]= j DC^[C,H+]. 

It is now just a matter of differentiating under the integral sign to show that the amplitude 
f ll65p does obey the Schwinger- Dyson equation (11611) . We have 

(166) fm[H+] = j DCf =ih* j = . 

A more rigor treatment of the (divergent) Gaussian integral (I165P implies fixing the gauge 
freedom ( 1164p by the BV method for reducible gauge-algebra generators [S], [6]. 
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